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integral equations of the second kind. The parameter estimation technique by 
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assumption imposed on the parameter space, the convergence analysis for the 
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here are applied to a boundary parameter estimation for 2-D elliptic systems. 
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1 . INTRODUCTION 


Recently, there Is growing interest in the problem of identification for 
distributed parameter systems (IDPS) both from the theoretical and numerical 
points of views (e.g. [4~8], [ 1 1 ] , [ 1 5 ] , [2 1 ] , [ 22] ) , The minimization of the 
output least square criterion (OLS) is one of the most popular methods for 
TOPS (see [16] and the references therein). The numerical determination of 
parameters by minimizing OLS involves some theoretical difficulties. Namely, 
numerical methods proposed must be equipped with convergence properties, such 
that the solution of discretized optimization problem implemented on the 
computer converges to the optimal solution of the original problem in infinite 
dimensions when we take the limit with respect to the number of dimensions. 
The compactness idea in the context of parameter estimation provides us a use- 
ful theoretical framework for the convergence and stability arguments in 
computer implementations of the discretized problems (see [4~8]). 

In this paper, an effort for parameter estimation is directed to the 
numerical method of parameter identification on the boundary, equipped with 
convergence properties for the optimal solution of discretized parameter esti- 
mation problem. Let 0 and 9G be the bounded domain and its boundary 
curve such that 


30 = {x = 5(t)|?(t) = (5 1 (t),5 ? (t)), 5 k (t) (k = 1,2) are a 
1-perlodic class function on [0,1] with |d£/dt| t 0} 

where { r } : = -[-r] and [r] stands for the greatest Integer _< r. 

Throughout this paper, we assume r 2. We consider the boundary integral 
equation of the second kind, 
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1 

<J>(s) - / 

0 


k(s,t;0)<Kt) |^|-|dt 


f(t) 


( 1 ) 


where 0 is the unknown parameter to be identified. A large number of 
classes of 2-D elliptic boundary value problems can be reduced to systems of 
Integral equations of the second kind as described in (1). Such integral 
equations are found in many applications to the field of engineering which 
appear in thermal diffusivity, viscous flow, electrostatics, acoustics, 
elasticity, etc. (See [10] and the references therein). 

Motivated by this fact, we consider the following parameter estimation 
problem. 


p 

(TP) Given the measurement data ^ at t^eTO, 1 ] (i = l,2,»«*,m), 
find 8 €® which minimizes the OLS 


•K8) = y I |4>(tP,0) - «K| 2 
i=l 1 1 


( 2 ) 


subject to (l). 

In Section 2, we discuss the existence property of an optimal solution for the 
problem (IP), In Section 3, the finite approximation technique for parameter 
estimation is proposed. Theoretical convergence proofs of the approximation 
method proposed are given. In Section 4, the method proposed in Section 3 is 
applied to the boundary parameter estimation of 2-D steady state thermal dif- 
fusivity, Finally, some numerical results will be demonstrated in Section 5, 
The notation used throughout this paper is standard and closely follows 
that explained in [17], For norms, we use !!•! where X is an appropri- 

A 
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ate Banach or Hilbert space. We denote by H r (Y) the periodic Sobolev space 

of arbitrary real order r on Y. The notation C^ c ^(Y) stands for the 

• • 

space of {t}-tlmes Holder continuously differentiable functions on Y. 


2. PRELIMINARIES 

The equation (l) is rewritten symbolically as 

f(I + K(0))<H(s) = f(s) for 0 < s < 1 (3) 

with the integral operator 


K(0) := - / k(s,t;0)(.)|^||dt. 

0 az 

Throughout this paper we assume that: 

(H — 1 ) 6€C^^(T;R^^) 

where T is a closed interval in [0,1], t is a real number such that 
t < r - -i- , and I is a natural number, respectively; 

(H-2) aH<J>ll < II (I + K(0 ) )4> I 

H r (3G) H r (3G) 


for each 0 and ViJieH r (3G); 
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(11-3) K(0 ) is a compact linear operator in H r (3G) 


for each 0; 


(H-4) 


1 (K(9 . ) - K(0,))*l < 0110 

1 * H r (3G) “ 



C* t} (T;R (A) 


1*1 r 

) H r (3G) 


for V<|> 6 H r (3G); 


(H-5 ) 


f €H r (3G). 


The admissible parameter set © in (IP) is then defined as follows: 


© is a compact subset of H r (T;R^) and its metric space (®,d) is com- 
pact with respect to the distance function 


d := p(x,y) = dx - yl 


C^CTjR^) 


for x,y£(S>* 


Theorem 1 : Under hypotheses (H-l) to (H-5), there exists at least one 

solution 0*€© of (IP). 


To prove Theorem 1, we need the preliminary results: 


Proposition 1: Let $ ^ ■ 4> (0 ^ ) and ^ * <|> (0 2 ) he the solutions of 
(2) corresponding to the parameters 0^ and 0^. Then, under (H-l) to 
(H-5), there exists a positive constant a such that, for V €®, 

V y 2 
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«♦, - ♦ ,« r < K p(6 ,0 )»fl 

H r (3G) a H (3G) 


( 4 ) 


Proof: From (3), It follows that 


(i + K(e 1 ))4 1 - * 2 ) - (K(9 ? ) - KCe 1 ))4> 2 - 


By taking the norm in H r (3G) and by virtue of (H-2) and (H-4), 


•♦l - <M r <1 r 

1 1 H (3G) 1 1 1 H (3G) 


From this fact and noting that 


‘V r <5- 1(1 + r " 5T ,fl r 

1 H r (3G) H (3G) H (3G) 


(5) 


we obtain the a priori estimate (4). 


Proposition 2: For (2), the following inequality holds: 


|j(0.) - J(e,)| <5^ i*(e.) - ♦<e,)i . 

1 1 “ Z 1 Z H 1 (3G) 


m 


X (114, (0 )| 2 + «<K0 )H 2 +4 \ |<|» | 2 ) 1/2 . 

1 H (3G) Z H (3G) i=l 1 


( 6 ) 


By virtue of Sobolev imbedding (see [17, p. 46]), we find that, for any 


i|>€ H l (3G) 


m 


l i , * I COO) <"V W) - 


By applying this, we obtain the inequality (6). 
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Proof of Theorem 1: From the compactness property, we may extract a sub- 

sequence {0^,} of a minimizing sequence {0 n > in © such that 0^, -»■ 0 
as n' > ®. From Propositions 1 and 2 and using (5), we have 


|j (0 ,) - j( 0 )| < p(e ,,e*)nf» 

n "a 2 n H r (9G) 


c{I nfB 2 +2 l l^l 2 } 1 ' 2 < c p (8 ,0*), 

“ H r (3G) i=l 1 n 


Hence we argue that 


J(e n ,) + J(e*) 


* 

inf J(0) as 0 ^ *■ 9 . 
0 € © 11 


The proof has been completed. 


3. PARAMETER ESTIMATION TECHNIQUE AND CONVERGENCE ANALYSIS 

In order to solve (IP) on the computer, we consider a discretized optimi- 
zation problem. Many numerical methods have been proposed for solving the 
Integral equation (e.g., see [ 1 ] ~ [ 3 ] , [20 ] ) • Here we adopt the collocation 
methods with spline functions* We select an increasing sequence of mesh 
points , 

tj - ^ - jh with j = 0,1,.. .,N; A N := {tj}" Q 


t = (.1 - j)h with j = 1,2,... ,N. 


and the nodal points 
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N 

We denote by S^(A ) the space of all 1-periodic (k-l)-times continuously 
differentiable splines of degree k with knot sequence A^. The spline 
collocation method for (3) is stated as follows [1], [20]: 

Find <(> N e S^(A^) such that the collocation equations 

[(I + K(0)H N ](t j ) * f(t ), j = 1,2,-. ,N (7) 


are satisfied. 


To formulate the discretized optimization problem, we approximate both the so- 
lution <J> and the admissible parameter 0 by a parabolic B-spline [9]. 
Let {^ 1=0 and ^ as ^ c elements with the knot sequence A N 

and A M (CT), respectively. Then we approximate <(i and 0 by 


M N+l MM 

♦ (t) = I wV(t) 

i=0 1 • 

£(M+2) 

0 M (t) = l xWt) 

i=l 


where 


(Bj, 0,0, -.,0) 


for i = 0, • • » ,M+1 


B j(t) 


(0 > ®i_M~2 »°» * * * » n) 


for i = M+2 , * * * , 2(M+1 ) 


(0,...,0,B^_ (ji _ 1)(m+2) ) for i « *M-M+£,«..,A(M+1). 
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The collocation equation (7) yields the linear equation for the unknown coef- 

N N N N 

ftclent vector w = col(w , • • • ,w^) , 

(C N + K N (9 M ))w N = f N (8) 


where C N , K N (0 M ) and are N-dimensional element matrices and vector, 

N 

respectively* Since <(> (t) claims a 1-periodic function on [0,1], we note 

N N N N M 

that = w 4t and w. Til = w,. We introduce the mapping l 11 of 

0 N N+l 1 

C*(T;R^^) into r(£( m+2 )) such that 1^0) is a compact subset of 

r (£(M+ 2)\ By using this, let us define the admissible parameter class (B^ 

for the discretized parameter estimation problem (AIP)^> M 



*(M+2) 

e M (t,* M ) = l A(t), 

i=i 

, M .,.M ,M .M N ^ T M_.« 

\ = coi(x 1 ,x 2 ,»** ,x £(M+2) ) e I ©1 


Thus we implement the parameter estimation numerically by solving the problem 


(AIP) n »M Find Q M (X^* ) € <3^ which minimizes the approximate OLS, 

A 


t N /c .Mv 1 v I.N. P .Mv . 1 2 

J (0 ) = -X I !♦ (t.,0 ) - *.| 

i=l 


(9) 


subject to (8). 


In the sequel, we assume that all hypotheses stated in Section 2 hold only for 
the case r = 2. Our next results play a fundamental role in establishing 


the convergence proof* 
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Lemma 1 ([20, p. 104, Theorem 3.5]): For each 9 eC^(T;R^^), the 

following error estimate holds: 


HO) - 4 , (e)| < Yhl<K9)l - 

H 1 (9G) ” H (3G) 


Lemma 2 (fl, p. 353]):' For any 0€H 2 (T;R (£) ) and e£(0,i), the 


re 


exists 0 **€S 2 (A^) such that 


-0 M » i _ ( 0 < 2 (o) 

C 1 (T;R U; ) H Z (T;R U; ) 


where W = max |t i+1 - t 1 | for eA M * 


Theorem 2: For each N and M, there exists an optimal solution of 

(AIP) N ’ M . 


To prove Theorem 2, we require 


Proposition 3: Let 0 £ = 0 M (X^)€(H^ (1 = 1,2). Then 


|( t» N (0.)l . (Yh + 1 )B f D , for 1 = 1,2, 

1 H 1 (3G) HOG) 


( 10 ) 


.N/„ v . N 


e 


^"(e,) - ♦ (e,)i . < (Yh + l ) n f n 9 n0, - 0 9 p , 

1 £ «* 1 /<\ J-I \ ^ /*\ 1 ^ r * 1 


H (3G) ot 


H (9G) 


c a (T;R (£) ). 


(ID 


Proof: Applying lemma 1, it follows that 
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ii<j> n ( 0 . )n , < H N (0.) - ♦(0.)i , + H<K0.)H , 

1 H (3G) H (3G) H (3G) 


( 12 ) 


< (yh + 1 ) (0 ) n „ 

1 HOG) 


i — (a) 

Since 0^£C (T;R V ), we note that, from (H-2), 


ii<j.(0.)ii , < i ii fn , 

H(3G) H(3G) 


(13) 


Combining (12) with (13), we obtain the inequality (10). By using the same 
procedure as in Proposition 2, it follows that 


1)^(0,) - <f> N (0 9 )n , < I- ne 

H (9G) 


" e2, cW 4 > 


ii <j> N ( q ) n . 

) 1 H(3G) 


(14) 


From the inequality (10) to (14), the inequality (11) can be obtained. 


Proof of Theorem 2: 

subsequence of {9 M (A^J)} 


Since I^tB> is compact, there exists a convergent 
such that 


eV.) . .V'") 

n 


N M 

where X 9 is the solution of 


inf J N (0 M U M )) = J N (0 M (X N » M )). 


e M €(B" 

Hence it suffices to show that, for each N and M, 0^(A^) J^(0^(X^)) 


is 


continuous. From Proposition 2, we have 
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lAeV/)) - AeV, 1 )) 


<"« ,Ae M (x?>> - ♦ M ce M cx5». . 

Z 1 HOG) 


,„.N, q M„M nn .2 ^ »aN/ a M,,M Ull 2 |2xl/2 

x (n4* (0 (x.))n . + ii<j> (e (x,))n , + 4 l |<i», | ) . 

1 H (3G) HOG) i-1 1 


Applying Proposition 3, we obtain 


|AAx”)) - Ae M <x“))| 


< ^ (Yh + i)n0 M (x”) - 0 M (x”)n . ( . 

~ 2 1 2 C 1 (T;R U) ) 


a 


x»fn {(Iiil) 2 »fn 2 + 2 l |<k|^} 

HOG) “ HOG) i=l 1 


m 


2 , 1/2 


(15) 


The statement of Theorem 2 directly follows from (15). 

The next theorem shows that the solution of (AIP)^>M converges to the optimal 
solution of (IP). 

MM ft 

Theorem 3: The sequence (0 (a 9 )} admits a convergent subsequence 

{0 (X 1 )} such that 


M M * 

J (9 lX J *)) •> J(0 ) = inf J(0) 

06 ® 


ns N ,M^ + ». 
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_ Nj ,M^ ^ «N. ,M^ 

Proof: Since X J €1 (B> , there always exists 0^ J £<fl) 

*^k M A ^1 *^k 

that X - 1 = I 0, J . By using this, it follows that 


M-V M k 

P0 M (X - 1 K ) - 0 B 




N ,M^ a N ,M^ a N ,Mj^ * 

< n0 M (x 1 k ) - e J k n . (SL) + H 0 , j - 0 i . . , . 

A C 1 (T;R U) ) A C 1 (T;R U; ) 


Applying Lemma 2, we have 


10 (X J ) - 0 n 

* 


C 1 (T;R^^ ) 


I A J *Sl 




< 61X10 J ''n 2 f v + ie j ’ ** - e*n 2 m • 

A H Z (T;R U; ) H Z (T;R U; ) 


Noting that ® Is a compact set in H^(T;R^^) and that if -► 0 as 


M L + the inequality (16) asserts there exists convergent subsequence 

K 


{0 


M (T N J’ Mk )> with 0 M (T Nj ’ Mk ) - 0* 


«. The remainder 


as N , * 

m w N.,M U * M N, ,M. * 

of this proof is to show J (0 (X 1 )) * J(0 ) as 0 (X J ) -*■ 0 , 


N , M 4- Applying Proposition 2, we derive 

I K 


|/(0 M (X 1 ,Mk )) - J(0*)| 


< |J N (0 M (X j> k )) - J N (0*)| + | J N (0 *) - j(0*)| 


<^H N (0 M (X Nj ’ Mk )) -*V). , 

“ H(3G) 


/.xN/ ft M rr N j ,M k xx .2 . ..N, ** .2 

x(ii(j> (e (X J ))u , + u 4> (0 ) u . + 


m 


H^SG) 


1 ■ 4 I l*ilV /2 


H ( 3 1» ) i=l 


such 


(16) 


half 
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y?m M * * 

+ ~ h n (0 ) - *(0 )ll , 

H (9G) 


x ( n <t> N (0 * >a 2 . + ii<t>(0*)n 2 1 + 4 l |* | 2 1/2 . 

H (3G) H (9G) i=l } 


By virtue of Lemma l and Proposition 3, we obtain 

|j N (9 M (r Nj,Mk )> - J(6*)| 

< Is (Th + I>I0 M (X J,Mk ) - 6*1 , , . 

“a^ C 1 (TjR U, > 


xjfn {(lAti) 2 n fn 2 + 2 l |* | 2 } 1/2 

H(3G) HOG) i=l 


. u /2wr f (yh+i) 2 + l ,„2 m 

+ h • — S- — II f II J { s II f II « ■ L I r J I 

H(9G) a H (9G) i=l 


♦* l U|V«. 


From the fact that h •> 0 as », we conclude that 


N M 

|J N (0 M (X j k )) - J(0*)| > 0 as Nj,^ 


The proof has been completed. 


4. APPLICATION TO PARAMETER ESTIMATION OF ELLIPTIC SYSTEMS 

We consider a 2-D steady state diffusion system 


Au = 0 


in G 


(17. a) 
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with the boundary condition 

?£ + c(u - g) - 0 on 3G, (17. b) 

d n 1 

+ 0.U = 0 on 3G 2 (17. c) 

where the domain G represents some geometrical shape of the material and one 
part of the boundary 3G^ means the surface of system structure. The given 
input g stands for the external excitation through the surface 3G^ and 
c is a known parameter on 3 G^. The problem considered here is to identify 
the parameter 0 on ® ne P oss *ble application in the problem above 

mentioned arises in material testing, in particular, in the area of 
nondestructive evaluation approach in space technology (see [13]). The 
parameter 0 on 3G^ represents the conductivity and/or heat transfer 
which characterize the structural flaws of unvisual section of the boundary 
3G. The decomposition of 3G into 3G^ and 8G 2 is taken as 
3G 1 = {£(t)|t € [0/t]} and 3G 2 = {?(t)|t € By using Green's 

formula, Eq. (17) is reduced to an integral equation on the boundary curve 
(see [10, Ch. 2] ,[12]) 

4>(s) - U(s))" 1 [f c(t)£n|c(s) - C(t)|<(.(t)|^||dt 
0 ac 

+ J l 0(tHn|5(s) - 5(t)|*(t)|g||dt 

C (18) 
+ / In * n l 5(s) " 5(t)|*(t) |^||dt] 

= -(X(s)) 1 / c(t)ln|g(s) - 5(t)|g(t) |^|.|dt 
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whore 

<j>(s) := u(5(s)) for 0 < s < 1. (19) 


In the above equation (18), A(s) is given by 


/tt for se[0,l]/{t^ =1 
A(s) := | 

for s€t* i = 1,2 ,•••,£ 


where Is the internal angle of the boundary at t*. By setting as 

k(s,t;0 ) = (X (s) ) 1 (c(t,0)£n|5(s) - 5(t)| + £n|5(s) - 5(t)|)|-^|-| 

and 

f(s) = -(X(s))" 1 / £n|5(s) - C(t)|i(t)|4||dt, 

0 

we have the representation (l) where 

^ ( c(t) for 0 < t < t 

c(t;0) := . “ 

0( t) for t <_ t <_ 1 

( c( t)g( t) for 0 < t < t 
g(t) := • 

,0 for t < t < 1 

respectively. Although each integral in (18) becomes singular, its limit is 

computable in terms of its principal value. Without loss of generality, we 


assume that 
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sup 

0<t,s<l 


|C(t) - 5(s)| < 1. 


Moreover, the following conditions are assumed: 


( A— 1 ) geH^Ojl) and supp(g)C3G 1 ; 

(A-2 ) c€H 2 (0,7) and 

0 < Oj _< c(t) ^ ^2 < " a.e. in (0,t); 
(A- 3 ) eeH 2 Ct,l) and 

0 < Hj < 9(t) _< n 2 < « 

0(7) = c(t) 0(1 ) = c(0) 

0'(t) = c'(t) 0'(1) - c'(0). 


Then Lt can be checked that the operator K(6) and the input f satisfy 
the hypotheses (H-2) to (H-5) in Section 2, Hence the results of Section 3 
can be applied to this problem. Since the unknown function 0 in (17. c) 

O 1 M 

takes Its value in H (0,t;R ), 0 (t) in Section 3 is simply rewritten by 


M+l 

e H (0 - I x75?<t>. 


i=0 


11 


In order to assert the assumption (A. 3), we require the following linear and 



-17- 


inequality constraints: 

xj - cOt) - ^ c'(T) 

x? - c( t) ♦ ^ c'(7) 

( 20 ) 

M h u 

X M = c(0) " 2~ C ' (0) 

X M+1 " c(0) + T~ C ' (0) 

0 < 0 1 £ x i £ s 2 < 00 for 1 “ 2,..»,M-1, 

where 

v, _ 1-t 

h M M * 

The problem (AIP)^*^ in this section is to solve 

min J N (A M ) = J N (A N,M ) (21) 

X M 

subject to 

(C N + K N (A M ))w N = f N (22) 


with the constraints (20), where 


x M = U*,a m 


r 


x M ) 
,A M+r * 
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The corresponding coefficient matrices C^, K^(X^) and vector f N in (22) 
are explicitly given by 

f B^(t 1 ) + Bjj +1 (t t ) for 1 = 1 , • • • » N ; j = 1 

[0%^ := < B^) for 1 - !,•••, N;j = 2,...,N-1 

V B JJ( tl ) + B Q<t 1 ) for 1 = 1 , — ,N; j = N 

zr—f T(t ,t;X M ){Bj(t) + bJ! (t)}dt 

X(t t ) 0 1 1 N+l 

for 1 = 1,»«» ,N; j = 1 

1 , 3 ~ M N 

£— J T(t ,f,x")B^(t)dt 

ACtj) 0 1 3 

for i - !,•••, N;j = 2,***,N-1 

- r - J T(t ,t;X M ){B”(t) + Bp(t)}dt 

x(t t ) 0 

for 1 = l,2,»»*,N;j = N 

f f N ] : = _ _L_ J S(t.,t)g(t)dt 

3 A ( t ^ ) 0 3 

for ] « 1 , 2| • • • ,N 



where 
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T(s,t;X M ) : = {c(t,X M )*n|5(s) -5<t)| 

+ | 77 £n|5(.s) -5(t)|) |^| 
S(s,t) : = ln|e(s) - 5(t)| |^|| 

and where 


c(t,X M ) : 


c(t) 


M+l 


for 0 < t < t 


l lS(t) for t < t < 1 

1=0 


For a fixed N and M, the necessary condition for the optimality of (AIP) N » M 
can be obtained. 

Theorem 4: Set the matrix H M and vector tj> as 


B X> + B N + l (t ?> 


for i *» 1, • • • ,m; i 3 1 


■""in * 


( B N N (tf) ♦ 


for i 3 3 


for i 3 1 9 • • • ,ra; j = N 


s- 


for j 3 1 , • • • ,N, 


* n M 

Then the necessary condition for X * to be optimal is characterized by 
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T v H tt N ’ H )-iv 1 k n (; n - m )]'„ n (x n ' m )(x" - i"' M > > o 

k=0 k 


wlie re 


[C N + K N (X N * M )]w N « f N 


[C N + K N (X N,M )]'v N - -(H N )'(H N w N (X N,M ) - i|») 


and where xj^ Is an any element of r(M+2) with the constraint (20) , 


Proof: As Is well-known, the necessary condition for X 

optimal is characterized by 


:n,m 


v/(X N ’ M ) . (X M - X N ’ M ) > 0 


for all X M with the constraint (20). We note that, from (9), 


NM NM N NNM 

[v x J (x )] k = (w x (n)'("T(H V (i n ) - ♦) 

k 


for k - 0 , 1 , • • • , M+ 1 


where 




Mv 

w (X ) 


for j = 1,2,--- ,N. 


The sensitivity equations for (20) with respect to X^ 


(23) 

(24) 

(25) 

to be 

(26) 

(27) 


are given by 
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tC N + kV>,w» - — [V. K N (X M )1w N (X M ) 
k *k 

for k = 0, 1 , • • • ,M+1 . 


(28) 


Hence, by Introducing the adjoint equation 

[C N + K N (A M )]'v N = -(H N )'(H N w N (A M ) - ^), (29) 

we can evaluate the gradient vector by 

[V x J N (X M )l k - (w N U M ))'[V a K N (A M )Kv N (A M ). (30) 

A k 

Setting as X M - X N,M in (30), the substitution of (30) into (26) yields 
the variational inequality (23). The proof has been completed. 

In the sequel, we consider numerical procedures for solving (AIP)^>^. 
From Theorem 4, we can compute the gradient of the cost function using (30). 
Hence many optimization techniques for the constrained problem are readily 
applicable to our problem (see [14], [18] and their references therein). We 
can represent the constraint (20) as the linear equality and inequality: 

A*A M - b” (31) 


where 


a”a m < b» 


(32) 
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'l 0 0 • • 

0 1 0 • • 

1 0 0 0 • • 
.0 0 0 • • 


0 0 0 

0 0 0 

0 1 0 

0 0 1 . 



c(t) + 


c(0) - 


c(0) + 


V ' (t> 


h M c'(t) 


ly^CO) 

2 

h M c '(°) 

2 



— 

1 



-*2 " 


-1 

1 

0 


-»1 

• 

n 

CM 

< 

-1 

• 

• 

• 1 
- 1 _ 

n 

CM 

H , 


0 


_- s i _ 


For the numerical results reported in this paper, we adopted the gradient 
projection method which is a particularly useful technique for the 
optimization problem with linear constraints as described in (31) and (32) 
(see [14] , [18] , [19] ). The iterative algorithm by the gradient projection 
method can be stated as follows: 

<<Numerical Algorithm>> 

Step 0: Fix the number of dimensions N and M for the problem (AIP)^»^. 

Set an initial value vector X^ ,M (0) satisfying (31) and (32) 


and i = 0 
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Step 1: 


Step 2: 


Step 3: 


N N M 

Compute the gradient vector V J (X * (i)) by 

A 


[V x J N (X N ’ M (i))] k - w N '[V x K N (X N,M (i))]'v N 

k 


for k ® 0, 1 , • • • ,M+1 


[C N + K N (X N ’ M (i))]w N = f N 


[C N + K N (X N,M (i))]'v N - -(H N )'(H N w N (X N,M (i)) - <|i)< 


If A^X * (i) < b 2 , set 


n N ’ H (i) - -v/(* K - M (i)) 


and proceed to Step 5; otherwise, proceed to Step 3. 


(33) 

(34) 

(35) 

(36) 


Compute the current direction by 


P M V x J N (X N ’ M (i)) 

|P M V x J N (X N,M (i))| 


where 


.M 



(AVS^Ap 

p p 


(37) 


and A includes the gradient of all currently active constraints 

M N M 

associated with matrix k^. If n * (i) * 0, proceed to Step 4; 
otherwise, proceed to Step 5. 



Step 4: 


Step 5: 


Step 6: 


-24- 


Compute satisfying 


J N « N - M <i> + ,<;V- M (i» - min . J H (X N - M (i) + Yn N,M (i)) 

" ln t€[0,Yl 


( 38 ) 


A N M 

where y is the largest step that may be taken from A 9 (i) 

N M /i\ * 

along n * (1) without violating any constraint. If Y v . / = Y» 

min 

then add the new constraint to the matrix A^ and proceed to Step 
5; otherwise, the new approximation to the solution is given by 


a N,M (1 + d = x N ’ M (i) + Y (i) n N,M (i) 

min 


(39) 


and proceed to Step 6. 


Compute the vector 6(A^’^) by 


6(A N » M (i)) = -(A M A M ') L A M V J N (A N,M (i)). 
P P PA 


(40) 


If all components of 6 are nonnegative, then set 


X N ’ M - A N ’ M (i) 


and terminate the computation; otherwise, delete the column of A 

N M 

corresponding to the smallest components of 6 (A *(!))• 


M 


If i < i , then replace i + 1 by i and return to Step 1; 
— max 

otherwise, print the statement "iteration over" and stop the 


computation 
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5. NUMERICAL EXPERIMENTS 

This section is devoted to a report of our efforts on computer implemen- 
tation of our techniques for the problem (AIP)^>M outlined in the previous 
section. In numerical experiments, we set the rectangular domain G with the 
boundary curve 5(t) given by 


'(2t, 0) 


for 0 < t < 7 - 
— 4 


€(t) - (5j(t), C 2 (t)) 


|( 2 * 2t “ 2 ^ 


| (“ 2t + * -j) 


for <. t < j 


for j < t < | 


k(0, -2t + 2) 


for ^ < t_< 1. 


Clearly, we have 


isii ■ < 4 r> 2+ (^) 2 ) i/2 - 2 £ ° r 0 < t < 1 


(42) 


_ 3 

and we decompose the boundary 8G into 3G^ and 3G 2 by t = -j- . 
The known function c(t) was preassigned as 


c(t) = 10 


for 0 < t < =r . 
— — 4 


(43) 


The test input was set as 


for 0 < t < ^ 


g(t) = 10 


(44) 
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The unknown function ©^(t) t * ie numer ^ ca ^ experiments was given by 


9 0 (t > = 0 O ( 4 " C) 


for t < 1 


(45) 


-800t + 1200t - 440 


3 31 

for l< * <& 


e 0 (t) 


-40t + f 


. 31 ... 17 

for 40 i ' < 20 


2 1237 

800t - 1400t + - 


£or §<'<? 


Solving the boundary element model (34) and (35), the computations of 
element matrix K N and vector f N must be carefully treated since, as 

indicated in Section 4, K N and f N involve the integrands with the weak 


singular kernel. In the numerical experiments, we computed these in the fol- 
lowing way. The ( i , j )-component of K N (i = l,2,***,N;j = 2,*«*,N-1) is 
rewritten as 

N t £ 

- £ [ ^ / c(t,X M )Jln|5(t 1 ) - 5(t)|Bj(t)dt 


'£-1 


(46) 


+ I / in|?(t ) - 5(t)|B^(t)dt] 

o=l t dn 1 J 

£-1 


For the first integral, we decompose it into 


N £ M 

I / c(t,X M )inU(t ) 

t £— 1 

£*i 


5(t)|B”(t)dt 
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+ / c(t,X M )£n|5(t ) - 5(t)|B I J(t)dt. 

t 1 J 

i-1 

The first Integrand becomes well-behaved and each integral can be approximated 
using an appropriate Gaussian quadrature formula. On the other hand, the 
second integral becomes singular. In this case, we further decompose it into 


'■1 

/ c(t ,X M )fcn| tj- t | B^( t)dt 

C i-1 


+ / 


iv M 

c(t,X Kn 
i-l 


|€Ct ± )-€(t) | 

\t t - t| 


Bj( t)dt. 


We can compute the first term analytically. The second integrand has a remov- 
able singularity at the point t^ since |-g|-| * 0. Therefore, the 
Gaussian quadrature can again be used to approximate this integral. Similar- 
ly, the second integral in (46), which is the so-called "double layer 
potential," has a removable singularity at point t^ (see [12] for more 
details). The element vector f^ can be computed in the same way. 

Our numerical experiments were carried out as follows: First, the ob- 
served data Y were generated by solving the boundary 

1 z m 

element model, i,e, , 




(i = ,m) 


M" M" W" 

[C N + K N (0)]* w = f 


with the true parameter 0 = 0^ for N' > N, Secondly, we implemented 
the numerical scheme proposed in Section 4. Experiments were carried out for 
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the dimension N ■ 40 and M = 10. Namely, we solved the problem 

(AIP)^ 9 ’*^. The initial guesses were given by 

[A 10 (0)]j ■ 10 for j ■ 1 , • •• , 10. 

In Step 1, the numerical Integration of and f^ can be accomplished by 

using the eight point Gauss-Legendre formula. In Step 2, parameters of the 
constraint vector b£ in (32) were chosen as 

« 2 S 2 = 12. 

In Step 4, the golden section search was used for a method of optimization 
along a line [14]. In Step 6, we set i max = 64. 

In the numerical example, we checked the sensitivity of our estimation algo- 
rithm with respect to the number of sensors by testing experiments for a 
different number of sensors. From several numerical findings, we suggest 

that, for a couple of fixed dimension (N,M) of (AIP)^ * , the proposed method 
requires at least M sensors, i.e., m J> M, in order to assert the 
effectiveness of the proposed method. To demonstrate this, we present 
numerical, results for the following three typical cases; 

(Case l ) m = 20 


rP = 21+9 

l 80 


for i = 1 , 2, • • • ,20 
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(Cnse 2) 


m = 10 


41+7 * 4 1 , 

"To" for 1 = 


33 

m 


for 1=6 


41+9 

80 


for 1 = 7, • • • , 10 


(Case 3) 


m = 4 


121-1 

80 


for 1 = 1,2 


121+1 

80 


for 1 = 3,4. 


The corresponding sensor locations are Illustrated in Fig. 1. Our results are 
given In Tables l, 2, and 3. For convenience of comparative discussions, the 


M 

estimated parameter function 0 (t,X n>11 ) and true function 0 n (t) are 




O' 


shown In Figs* 2, 3, and 4, respectively. 


6. CONCLUDING REMARKS 

The computational costs in solving the identification of distributed 
parameter systems often tend to be quite expensive. For some class of 
boundary value problems (BVP), computational savings can be achieved in the 
use of the boundary integral equation method (BIE). Since the BIE method 
makes it possible to replace the BVP by certain integral equations, the 
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application of the BTE method to the related identification problems has the 
effect of reducing the dimension of the problems. A couple of effective 
estimation algorithms by B1E methods have been proposed for identifying 
boundary parameters without the theoretical convergence proofs [211, [22]. In 
this paper, we developed the numerical method with emphasis on a theoretical 
framework for the boundary parameter estimation technique. To assert the 
convergence property, we claim the regularity of the solution for an integral 
equation model and also require compactness of the set of parameter functions 
to be identified. A spline collocation method of even degree was used to 
obtain a spline-based parameter space. Following the compactness ideas for 
parameter estimation in [4]~[8], a convergence property for the discretized 
solutions of (AIP)N>M was shown in Section 4. The efficiency of the theo- 
retical idea was demonstrated by the numerical experiments. Our approach can 
be readily extended to treat optimal shape design, optimal shape control, etc. 
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Table 1. Estimated Value of 


X N ’ M 


in Case !• 



initial guess 

iteration 8 

iteration 16 

iteration 32 

\ N,M (iteration 64) 


10.000 

8.7072 

8.7377 

8.8241 

8.8812 

\ 

10.000 

7.7119 

7.7159 

7.7352 

7.7494 

A 4 

10.000 

6.9597 

6.9432 

6.9102 

6.8905 

A?’" 

10.000 

6.5553 

6.5277 

6.4654 

6.4263 



10.000 

6.5553 

6.5267 

6.4655 

6.4262 

\ 

A r 

10.000 

6.9539 

6.9406 

6.9106 

6.8902 

^8 

10.000 

7.7028 

7.7117 

7.7360 

7.7491 

\ yv,M 

10.000 

8.6924 

8.7319 

8.8250 

8.8807 

Value of 

Cost 

Function 

6.1742 

XlO -8 

— 

'■gf ■ 

■ 

1 

7.6534 

xlO -8 

5.2906 

xlO -8 



t 


Fig. 2. True and estimated function in Case 1 
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Table 2. Estimated Value of 


in Case 2. 



initial guess 

iteration 8 

iteration 16 

iteration 32 

\N,M (iteration 64) 

\ w,AI 
* 2 

10.000 

7.8689 

8.0285 

8.1654 

8.3297 

A 3 

10.000 

7.5169 

7.5837 

7.6076 

7.6438 

\ nM 

A 4 

10.000 

7.2489 

7.2387 

7.1752 

7.1108 

\ N * M 

A h 

10.000 

7.1080 

7.0479 

6.9372 

6.8164 

\N,M ~ 

10.000 

7.1215 

7.0452 

6.9370 

6.8159 


10.000 

7.2887 

7.2311 

7.1749 

7.1094 

A s 

10.000 

7.5817 

7.5716 

7.6076 

7.6421 

A 9 

10.000 

7.9562 

8.0133 

8.1664 

8.3289 

Value of 



5.4085 

4.1514 

2.7980 

Cost 

Function 

HI 


xlO -7 

xlO** 7 

XlO -7 



t 


Fig. 3. True and estimated function in Case 2 
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Table 3. Estimated Value of 


X N » M in Case 3. 



initial guess 

iteration 8 

iteration 16 

iteration 32 

X N,M (iteration 64) 

-rnja 

A i 

10.000 

7.2629 

7.7044 

7.7634 

7.9721 

% yV'M 
A 3 

* 10.000 

7.4443 

7.5174 

7.5326 

7.5681 

-jim 

10.000 

7.5814 

7.3820 

7.3639 

7.2961 

A?’ M 

10.000 

7.6516 

7.3108 

7.2726 

7.1076 

\N,M 

10.000 

7.6511 

7.3110 

7.2722 

7.1073 

A 7 

10.000 

7.5800 

7.3853 

7.3626 

7.2684 

~ 

A 8 

10.000 

7.4419 

7.5227 

7.5307 

7.5673 

A d 

10.000 

7.2596 

7.7119 

7.7611 

7.9714 

Value of 

Cost 

Function 

2.3216 

xlO" 8 

4.8032 

xlO" 7 

2.8174 

xlO -7 

2.6762 

xlO -7 

1.8701 

xlO” 7 



t 


Fig. A. True and estimated function in Case 3 
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